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Abstract 

An analytical model was developped to describe the current induced DW dynamics of a Bloch 
DW in the presence of an external transverse magnetic field. The model takes into account the 
DW deformation and the magnetization tilting in the domain. The model is compared to the 
results of micromagnetic simulation and an excellent agreement is obtained. In the steady state 
regime, the model shows that the domain tilting does not change the DW mobility. An external 
or current induced transverse magnetic field such as the Oersted or Rashba field can prevent the 
Walker breakdown leading to a higher domain wall velocity. 
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INTRODUCTION 



The possibility to manipulate a domain wall (DW) using a spin polarized current has 
opened a new path to write the information in magnetic memory devices Several in- 
novative memory devices based on this effect have been recently proposed. One can cite 
three-terminal MRAM schemes where the different writing and reading path potentially 
solves the reliability issue observed in standard STT MRAM {2, 3]. The magnetic racetrack 
where the information is stored in tens of DWs simultaneously moved by current pulses is 
also a potential disruptive technology for the mass storage market which offers a fast, cheap 
and solid state alternative to hard disk drives [4] . The key factors for such memories are the 
write(/read) time and the power consumption of the device. An important effort of research 
has been carried out these last years to find out the ideal material with fast DW motion 
(> 100 m/s) and low critical current density (~ 10 11 A/m 2 ). While most studies initially 
focused on the soft materials such as permalloy, the attention has now shifted to out-of-plane 
magnetized materials with strong perpendicular anisotropy where spin transfer was shown 
to be more efficient. These materials also have the advantage of narrow DWs (~ 5 nm) 
meaning a potentially very high areal density. Besides the material issues, the nature of the 
DW dynamics plays also a major role. Indeed, two dynamical regime are expected depend- 
ing on the current density |5| : at low current density, a steady state regime with a stable 
DW structure and at high current density, a turbulent regime where the DW structure con- 
tinuously oscillates. The limit between these two regimes is the so-called Walker breakdown 
similar to the one observed when applying an easy-axis magnetic field. In large perpendicular 
anisotropy materials where high non-adiabatic effects have been identified one expects a 
much higher DW mobility in the steady state regime compared to the precessional regime. 
However, a current density higher than the Walker current is generally needed for fast DW 
motion so the low mobility turbulent regime is reached. Recently, several schemes 6|-|9| have 
been proposed to prevent the occurrence of this low mobility regime using in-plane magnetic 
field transverse to the wire where the DW propagates: the transverse field stabilizes the DW 
structure and thus extends the high mobility regime toward higher current density. This 
transverse magnetic field can be externally applied (9] or more advantageously induced by 
the current injection using the current Oersted field in multilayer geometry jf], Q, [lO] or 
the Rashba field in inversion asymetry multilayer structure jll]. Whereas the effect of a 
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transverse magnetic field on the field driven DW dynamics has been studied both experi- 
{9! I12I-2II and theoretically 22-25], the effect on the current driven dynamics has 



mentally la 112 

been little addressed. Of particular interest is the case of out-of-plane magnetized materials 
where the DW dynamics can generally be well described by a ID analytical model due to 
the simple Bloch DW structure. 

Here we present an analytical model to describe the current induced DW dynamics of 
a Bloch DW in the presence of an external transverse magnetic field. The model is based 
on a modified Bloch magnetization profile that takes into account the DW deformation 



and the magnetization tilting in the domain 



251 ] . The model gives a simple analytical 



expression of the current induced DW velocity in the presence of a transverse magnetic field 
in the steady state regime and an excellent agreement with micromagnetic simulations is 
obtained. The model shows that the transverse magnetic field leads to an increase of the 
Walker current allowing high mobility DW motion even at high current densities. 

The one-dimensional model initially proposed by Slonczewski [26[ and later extended to 
include the effect of spin torque jsl, Q, provides a convenient and intuitive way to describe 
the DW dynamics in the presence of current. This model is based on the assumption that 
the DW keeps its static structure during its motion determined by the exchange and the 
anisotropy whereas the forces inducing the dynamics are first order correction affecting 
only the DW position q, the DW angle ip and the DW width A. This model was shown 
to provide a good qualitative description of the DW dynamics of the Neel wall found in 
in-plane magnetized wires and a nearly quantitative description for Bloch wall 28 ]. 

The dynamical ID equations are based on the Landau-Lifschitz equation where the cur- 
rent induced spin torques have been added : 

dm 5E dm dm dm 

^ = - 7 5M Xm + " mX ^-^ + /3umX ^ W 

Here 7 = g\/j,B\/h is the gyromagnetic ratio and E the free energy density. The third term 
is the adiabatic spin-transfer torque where u = JPgfiB/2eM s , \x& the Bohr magneton, M s 
the saturation magnetization, J the current density. The fourth term is the non-adiabatic 
torque [s| described by the parameter (3. 

In the following, two cases are distinguished depending on the amplitude of the external 
transverse magnetic field H t compared to the anisotropy field H an = 2K an / fi M s . 
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Low transverse magnetic field 



For H t -C H an , the DW structure and the magnetization in the domain are little affected 
by the transverse field and one can assume a standard Bloch profile. The polar and azimuthal 
angles 9 and (p are then assumed as 9 = 2 arctan(exp{[x — q(t)}/ A(ip)}) and f = ip(t) with 
ip(t) constant. Here A(ip) = \J A/k is the DW width with k = (K an + K sin 2 ip) 1 / 2 where 
K an = K® n — fioM 2 with K® n the easy axis magnetic anisotropy, A the exchange constant, 
K = fioM s Hk/2 with Hk the DW demagnetizing field. The integration of Eq. [I] over this 
DW profile leads to the following ID equations: 



— {q u) +ip 

A a 

q — u — aAtp 



7 da 
~2M s ~dq 
7 da 



2M s dijj 

For H t <C H an , the DW energy per unit surface can be written as 

a = aVAk + H k M s A sin 2 ip - 2HM s q - irAM s H t cos(^ - ip t ) 



(2) 
(3) 

(4) 



with H the easy external magnetic field and ipt the in-plane direction of the transverse 
magnetic field 3l|. Eq. H]then leads to the following ID equations : 



^(q-P u )+iP = jH (5) 

A ; l^H k . ^AH t . 
q — it — uAip = — - — sin Zip H sm(^ — ipt) (o) 

The transverse magnetic field thus enters the equation as an additional torque on the DW 
proportional to sin(ip — ipt)- 



Large transverse magnetic field 

For larger H t , the magnetization tilting in the domain as well as the DW structure 
deformation have to be taken into account. This problem was already tackled in the case of 



field driven DW dynamics by Hubert 22j, [29]. To account for the tilting of the magnetization 
in the domain and DW deformation, more convenient polar and azimuthal angles 9 and (p 
are defined (see Fig. [[]). Their correspondent axes are rotated by an angle \ that defines 



the tilt of the DW plane due to DW dynamics. The application of the transverse magnetic 
field along a direction ipt leads to a tilt of the magnetization in the domain along 9 = 9q and 
if = ±<p with 



cos 9 Q = h sm(ip t — x) 
sin 9 cos ifo = h cos(ipt — x) 



(7) 
(8) 



where h = H t /H an . In between, the angle if rotates as |22|, |25j : 



cos <f — cos <f 



sin 2 (fo 



cos f + cosh[k (x)(x - q(t))] 



k {x) = sinyW {K an + K x)/A 



(9) 

(10) 
(11) 



The application of the transverse magnetic field leads to a widening of the DW (see 
Fig. DJc)) with a DW width scaling with h as 22j : 



A t (h) 



[y/l-h 2 sm 2 (4> t - x) - h cos(^i - x)\ \A + K d/ K an sin 2 x 



(12) 



To describe the current induced DW dynamics for such a DW profile, a Lagrangian approach 
is considered. Indeed, the LLG equation can be derived by applying the Lagrange-Rayleigh 
3 27, 30] 

f)T, ft fti, ft Pit, ftp 

(13) 



equation 



dL d dL d dL dF _ 
dX ~ didX ~ dxdX' + dX ~ ° 



with X J =dX/dx and using the Lagrangian Lt : 



L t = E + —(cos 9- cos^o) - (uM s /-f)fd( cos 9) /dx 

7 



(14) 



where E the micromagnetic energy density and 9 and f the polar and azimuthal angles de- 
fined as m = (cos 9, sin 9 cos if, sin 9 sin if). The angles 9q ans ifo describes the magnetization 
tilting in the domain with cos^o = hsimp t and sin^ocos^o = hcosipt- The Lagrangian is 
chosen so that it cancels out in the domain. The effect of the damping and the non-adiabatic 
torque is included in the dissipative function F [3^] : 

(15) 
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, d ^ f3u d 
v dt a dx 
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The integration of L t and the dissipation function F over the DW profile of Eq. [9] leads 
to |22| : 



LDWt 



a 



{2M s /j)(q-u)f k 



aM s E G 
7 4A 



ft 



a 



k 2 
h 



( dl/h 



„i a l 



V dx 



(2 ~ 2 
(7T -if - 



2<p 2 



1 — (/? cotan<£> c 



where 



fk 



•ri) 



cos( 



cos 9 )dip 



Y?0 



(16) 
(17) 

(18) 



E G = A^AK^J(1 - h*sm 2 (ip t - X ))(l + ^ /^anSin 2 X ) 



\J\ — h 2 — h cos(ipt — x) arctan 



(19) 



hcos(i/) t ~ X). 

Note that the adiabatic torque leads to the operation q — > q — u in the Lagrangian whereas 
the non-adiabatic torque leads to the operation q — > q — ftu/a in the dissipation function. 
The Lagrange-Rayleigh equation with X = x or q leads to the following ID equations : 



Eg. 



ft 



df k . 



a-^{q--u) + ^x = lHVl^ (20) 
4A a ox 



dfk,. , 



atE, 



G 



AA 



1 /3l/fc \ 2 l 



;.2 
a 



/ 2 ~ 2 
(7T -yj 



1 — (/?ocotan<^o 
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G 



2M S d X 



(21) 



From Eq. [20], one can easily derive the DW velocity in the steady state (x — 0): 

aE G 



(22) 



to be compare to the steady state DW velocity in the absence of transverse magnetic field 
obtained from Eq. [5] : 



7A 

q = ftu/a H H 

a 



(23) 



An important result is thus that the DW mobility induced by the current (ftu/a) is not 
affected by the transverse magnetic field in the steady state regime. On the contrary, the 
DW mobility induced by the easy axis field H is increased when applying the transverse 
magnetic field due to the large decrease in the DW energy Eq when H t increases. 
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FIG. 1: (a) Schematics of the DW profile, (b) Coordinate system related to 9, (p and \- ( c ) 
DW width as a function of the transverse field amplitude applied along the y direction (ipt = 0). 
The following magnetic and transport parameters are used corresponding to a Pt (3 nm)/Co (0.6 
nm)/A10x mutltilayer [g| : M s = 1.09 x 10 6 A/m,A = 1.6 x 10~ n J/m 2 , K° an = 1.25 x 10 6 J/m 3 , 
H k = 15 mT. 




Current density (1 12 A/m 2 ) Transverse magnetic field (T) 



FIG. 2: (a) DW velocity as a function of the current density for different transverse magnetic field 
applied (tpt = 0) predicted by the ID model (plain lines) and micromagnetic simulation (square 
plots), (b) Walker current as a function of the transverse magnetic field predicted by the ID model 
(square) and micromagnetic simulation (line). The following magnetic and transport parameters 
are used: M s = 1.09 x 10 6 A/m,A = 1.6 x 10~ n J/m 2 , K° an = 1.25 x 10 6 J/m 3 , H k = 15 mT, 
a = 0.5, p = 1.7. 

To obtain the DW velocity above the Walker breakdown, numerical resolutions of Eq. [20] 
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and Eq. [21] have been carried out (Fig. [5(a), plain line). The main effect of the application 
of the transverse magnetic field (ipt = 0) is an increase of the Walker current: the transverse 
field reduces the DW tilt angle and thus shifts the Walker current toward high current. Below 
the Walker current, the velocity follows v = (3u/a as expected, whereas in the precession 
regime well above the Walker breakdown, the velocity gets close to the zero transverse field 
velocity < v >= (1 + aj3) /(l + a 2 )u. The increase of the Walker current with the transverse 
magnetic field is shown on Fig[2^b). 
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FIG. 3: (a) DW velocity as a function of the injected current density at zero external transverse 
field and in the presence of a current induced transverse field (0.1 T for 10 12 A/m 2 ) along ^>t=0. 
The other parameters are the same as Fig. [2j 

To test the validity of the model, micromagnetic simulations have been carried out using 
the micromagnetic solver DWmicro3D based on the finite differences approximation for a 
100 nm wide nanowire and 0.6 nm thick Co layer 28| . A high perpendicular anisotropy is 
assumed leading to an out-of-plane magnetization and narrow Bloch DWs. The calculated 
velocity as a function of the current density and the Walker current are plotted in Fig. [2] 
(square symbols). One can see that an excellent agreement is obtained between the ID 
model and the micromagnetic simulations. Such a good agreement is generally not observed 
for in-plane magnetized Neel DW and can be accounted for to the rigid Bloch DW structure 
in this narrow wire. The profile of the magnetization transverse to the wire shows little 
dependence of the magnetization on the direction transverse to the wire. The DW thus 
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behave effectively as a rigid one-dimensional object. 

The model can also be used to simulate the effect of a current induced transverse magnetic 
field, such as an Oersted field or a Rashba field. Fig. [3] shows the velocity as a function of 
the current density for a transverse magnetic field of 0.1 T/10 12 A/m 2 . As expected, the 
current induced transverse field stabilizes the DW structure and thus prevents the Walker 
breakdown. The DW stays in the higher mobility steady state regime and a higher velocity 
is obtained. 

In conclusion, a one-dimensional model of the current induced domain wall dynamics in 
the presence of a transverse magnetic field that takes into account the DW deformation and 
magnetization tilting in the domains was developed. In the steady state regime, the model 
shows that the domain tilting does not change the DW mobility. The transverse magnetic 
field leads to an increase of the Walker current and thus the high mobility regime is extended 
toward high current densities. A current induced transverse magnetic field such as Oersted 
or Rashba can prevent the Walker breakdown leading to a higher domain wall velocity. This 
work was financially supported by the ANR 11-B510-008 ESPERADO project. 



Calculation of the kinetic energy 



fk is proportionnel to the kinetic energy of the DW E^ n 22J with : 



2M s g , s 

Ekin = -h A.24 
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It can be expressed in the (9, ip) coordinate system as: 



fk— (cos 9 — hsimjj t )(cosx + g(9, ^p) sin x)d(f (A.25) 



o 



where 



g(9,ip) = (cos<^ sin cosx — cos 6*o sin x) cos6'/(l — cos ^ ) (A. 26) 

cos# = cos <p sin 9 sin x + cos 9 cos x (A. 27) 

cos 0q = hsm(ip-x) (A.28) 

tan^o = T-^rf yT (A.29) 

hcos(i) - x)) 

For numerical calculation, an analytical expression of fk was obtained using the commercial 
software Mathematica. 
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